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ABSTRACT 

We study the late evolution of solar metallicity stars in the transition region between white dwarf 
formation and core collapse. This includes the super-asymptotic giant branch (super-AGB, SAGB) 
stars, which have massive enough cores to ignite carbon burning and form an oxygen-neon (ONe) 
core. The most massive SAGB stars have cores that may grow to the Chandrasekhar mass because 
of continued shell-burning. Their cores collapse, triggering a so called electron capture supernovae 
(ECSN). 

From stellar evolution models we find that the initial mass range for SAGB evolution is 7. 5. ..9. 25 Mq. 
We perform calculations with three different stellar evolution codes to investigate the sensitivity of 
this mass range to some of the uncertainties in current stellar models. The mass range significantly 
depends on the treatment of semiconvective mixing and convective overshooting. To consider the 
effect of a large number of thermal pulses, as expected in SAGB stars, we construct synthetic SAGB 
models that include a semi-analytical treatment of dredge-up, hot-bottom burning, and thermal pulse 
properties. To calibrate the time-dependent synthetic model, we have calibrated a number of SAGB 
stellar evolution models. This synthetic model enables us to compute the evolution of the main 
properties of SAGB stars from the onset of thermal pulses until the core reaches the Chandrasekhar 
mass or is uncovered by the stellar wind. Thereby, we determine the stellar initial mass ranges that 
produce ONe-white dwarfs and electron-capture supernovae. The latter is found to be 9.0. . .9.25 
for our fiducial model, implying that electron-capture supernovae would constitute about 4% of all 
supernovae in the local universe. Our synthetic approach allows us to explore the uncertainty of this 
number imposed by uncertainties in the third dredge-up efficiency and ABC mass loss rate. We find 
that for both processes, the most optimistic approach leads to about a doubling of the number of 
electron-capture supernovae, which provides a firm upper limit to their contribution to all supernovae 
of -20%. 

Subject headings: stars: AGB and post-AGB — stars: evolution — supernovae: general — stars: 
neutron 



1. INTRODUCTION 

It is well known that, for a given initial chemical com- 
position, it is the initial stellar mass which essentially 
determines the final fate of a star: lower masses pro- 
duce white dwarfs, higher masses neutron stars and su- 
pernovae. The late evolution phases of stars in the tran- 
sition region between white dwarfs and neutron stars is 
numerically difficult to model, and the relevant physics 
is not yet fully understood. This mass range is therefore 
often omitted in stellar evolution calculations. This is 
unsatisfactory because the uncertain initial mass region 
for this evolution is 7 to 12 Mq, implying that as much as 
half of all supernovae may originate from this transition 
region. 

Of particular interest in this context is the evolu- 
tion of so called super-asymptotic giant branch (SAGB) 



stars, which ignite carbon non-explosi vely, but also un- 
dergo thermal pulses ( Ritossa ct a l.1 Il996t llben et al.l 
1997t iGarcia-Berro et al.lll997l : iRitossa et al.l ll999': 'Sied 



20061) . These stars may end their lives either as massive 



ONe-white dwarfs ()Nomotolll984[ ) , or as electron capture 
supernovae (ECSNe), where the core collapse i s triggered 
by electron captures befo re Ne ignition (Wh eeler et al.l 
Il998t IWanaio et al] 12003 V Stars of larger initial mass 
ignite hydrostatic neon burning, form an iron core, and 
lead to classical core collapse supernovae (CCSNe). 

The upper mass limit of SAGB stars is affected by 
the second dredge-up, which may occur after core He- 
exhaustion, and which drastically reduces the mass of 
the helium core. At this point and throughout the follow- 
ing early SAGB phase carbon burning transforms the CO 
core into an ONe core (lNomotolll987t (Ritossa et al. !ll996t 
Iben et al.lll997l : iGarcia-Berro et al.i ri997: Ritossa 'etaLl 



Box 80000, NL-3508 TA 
a.j.t.poelarends@astro.uu.nl, 



MS B227, Los 
NM 87545, USA, 



^ Astronomical Institute, P.O. 
Utrecht, The Netherlands, 
n.langer® astro. uu.nl 

^ Theoretical Astrophysics Group, T-6, 
Alamos National Laboratory, Los Alamos, 
aheger@lanl.gov 

^ Keele Astrophysics Group, School of Physical and Geo- 
graphical Sciences, Keele University, Staffordshire, UK, fher- 
wig@astro.keele.ac.uk 

Department of Astronomy and Astrophysics, University of 
California, Santa Cruz, CA 95064, USA 

^ Department of Astronomy and Astrophysics, University of 
Chicago, 5640 S. EUis Avenue, Chicago, IL 60637, USA 



T999ir Since the temperature is not high enough to ig- 
nite Ne, the core cools, electron degeneracy in the core 
increases, and the structure of the SAGB star outside of 
the ONe-core is yery similar to th e most massive CO- 
core AGB stars (iFrost et al.lll998[ for a general review 
of AGB evolution see llben fc Renzhfil[l98aiHabing|[l99fil 
and Herwig 200^). The degenerate core is surrounded 
by the He- and H-shcU sources, which eventually pro- 
duce therma l pulses due to th e instability of the helium 
shefi source (jYoon et al.ll2004( l. 
In this situation, the mass of the H-free core contin- 
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ues to grow. If the core mass reaches the Chandrasekhar 
mass of 1.375 Mq, the core will collapse triggered by elec- 
tron captur es on ^'^Mg arid ^°Ne , and the star will become 
an ECSN (iMivaii et al.l [l980l: iMivaii fc Nomotol HqSTI: 
iHashimoto et al. logs'). Rece nt stu d ies bylRit ossa et al.' 
(119961): iGarcia-Berro et a.1.. dlQQTh: llben et al.. (,1997.) : 
iRitossaet al.l (|1999[ ) and lSiessI (|2006[ ) have shown that 
the mass fraction of ^^Mg in the ONe-core is smaller than 
previously thought, which diminishes the r o le of e lectron 
captures on ^**Mg. While [Gutierrez et al.l (|2005f) found 
that unburnt carbon in the degenerate ONe core could 
trigger an explosion at densites of ^ lO^g cm~^, we disre- 
gard this possibility furtheron since its observational im- 
plications are not worked out, and therefore this scenario 
can not yet be confronted with supernova observations. 

The initial mass range for core collapse after SAGB 
evolution depends on the effective core growth and mass 
loss of the SAGB star. Larger mass loss rates lead to 
a shorter duration of the SAGB phase. For very high 
SAGB mass loss, there is no time for any significant core 
growth, and the initial mass range for ECSNe will be 
very small. On the other hand, the core growth rate 
in SAGB stars depends on the the hydrogen shell burn- 
ing and thus on two crucial factors, hot-bottom bu rning 
(jSackman n fc Boothrovdlll991l:IVentura et al.ll2005D . and 
the efficiency of the third dredge-up. 

Previous studies of SAGB stars hav e concentrated 
on t he evolution of the stellar cores (|Nomotol 11984 , 
Il987| ). According to these models, stars with helium 
cores between 2.0 and 2.5 form ONe cores and ex- 
plode as ECSN, leaving a neu t ron st a r less massive than 
l.S Mr^. iRitossa et alJ (Il99l Il999l). llben et all (|l997h 
and IGarcia-Berro et aTl ()1997f ) studied the evolution of 
complete SAGB stellar models in detail. They describe 
SAGB thermal pulses, and an outward mixing event 
which they called dredge-out, in which the convective 
envelope connects t o a convecti on zone on top of the he- 
lium burning layer. ISiessI (|2006D . who studied the effects 
of the carbon flame and of thermonuclear reactions on 
the structure of the ONe core, finds similar results. 

Thermal pulses in AGB evolutionary models require 
high numerical resolution, both in time and space. The 
interpulse period decreases with increasing core mass to 
eventually only a few years for the most massive AGB 
star. In order to follow the evolution of SAGB stars with 
very high core masses, orders of magnitude more thermal 
pulses (up to 10000s) have to be computed compared to 
low-mass AGB stars, which experience only tens of ther- 
mal pulses. For this reason, no detailed stellar evolution 
calc ulations throu gh the entire super-TP- AGB phase ex- 
ist. IRitossa et al.l (|1999|) attempted to characterize stars 
that would end as ECSN. Based on the assumption of 
a constant SAGB mass loss rate of 10"'' Mq/ yr, they 
speculated that out of their set of four calculated models 
(9, 10, 10.5 and 11 Mg) only the 11 Mq model would ex- 
plode as an ECSN. The other models would lose all their 
envelope before the core has grown eno ugh, and their fi- 
nal fate would be an ONe white dwarf. lEldridge fc Tout) 
()2004br ) determine a minimum mass for supernova explo- 
sion around 7M0 (with overshooting), or around 9M0 
(without overshooting), again without being able to cal- 
culate the stellar evolution models through the final 
phases. 

Models of SAGB evolution suffer from two dominant 



sources of uncertainty: mass loss and the efficiency of the 
third dredge-up. To explore these uncertainties would 
require to compute several model grids, which is not fea- 
sible at this time. We therefore take a different approach 
and use the fact that TP- AGB stars, after a brief tran- 
sition phase, reach a quasi-steady state in which the im- 
portant structural quantities evolve in a simple and pre- 
dictable way as function of time. This approach of syn- 
thetic AGB modeling has already been s uccessfully used 
for low-mass a nd massive AGB stars (iRenzini fc Volil 
[19811 llben fc^enzini 198llMarigo et al.lll996D . 

In the following, we first describe the detailed stellar 
evolution models (Sect. [2]) and identify the initial mass 
range for SAGB stars by calculating the pre-AGB evo- 
lution phase up to the end of the second dredge-up and 
possibly C-ignition, using three different stellar evolu- 
tion codes (Sect. [3]). Next, we describe our SAGB steUar 
evolution models (Sect.|4]), and our synthetic SAGB evo- 
lution model (Sect. [5]). We present our results in Sect.|6| 
and concluding remarks in Sect. [71 

2. NUMERICAL METHODS 

We use three different stellar evolution codes to cal- 
culate the evolution of solar metallicity stars up to 
the end of the second dred ge-up, or to Ne-ignition. 
We u sed the codes ST ERN |Langed |1998|; iHeg er et a l.| 
l2000l). EVOL (lBlockedll9 95: Hcrwig 2000) and KEPLER 
(jWeaver et al.lll978l : iHeg cr ct al. 200 Q). AU th ree codes 
use the OPAL opacities (Iglcsias & R ogersll996[ ) , and are 
equipped with up-to-date input physics, including a nu- 
clear network with all relevant thermonuclear reactions. 

For our investigation, the most relevant difference be- 
tween the codes concerns the treatment of convective and 
semiconvective mixing. As we will see, these affect the 
He-core mass after central He-burning, and thereby the 
final fate of the stellar model. STERN and KEPLER use 
the Ledoux-criterion to determine convective instability, 
and take semiconvection into account. In KEPLER the 
treatment of semiconvection leads to rather fast mixing. 
Specifically, it is approximated as a diffusive process with 
a diffusion coefficient that is 10 % of the radiative diffu- 
sion coefficient. Typically, this leads to a lOOOx shorter 
mixing time scale as for the default value oflLange r et al.l 
(|1983| ) as used in STERN (asom = 0.01). No modification 
to the temperature gradient is assumed, i.e., the radia- 
tive temperature gradient is used. Additionally, in KE- 
PLER convection zones are extended by one extra grid 
point where fast mixing is assumed, to mimic convec- 
tive overshooting. In the EVOL code, convective bound- 
aries are determined by the Schwarzschild criterion, and 
semiconvection is not treated as a separate mixing pro- 
cess. Mixing beyond convective boundaries is treated by 
adopting an exponentially decayi ng diffusion coefficient 
()Herwig et al.lll997l : lHerwiell2000f ). Such mixing may b e 
induced by convective overshooting (Frevtag et al!'l99(}), 
or internal gravity waves ( Dcnisscnkov & Tout 2003) , or 
a combination of both (Yo ung et al.ll2005D . For the pre- 
AGB evolution, the overshoot parameter in EV OL has 
been set to / = 0.016, which was shown by iHerwigl 
([2000') to reproduce the observed main sequence width in 
the HR diagram of young open clusters. Effectively, the 
strength of mixing in KEPLER lies in between that of 
STERN (slow semiconvective mixing) and that of EVOL 
(Schwarzschild criterion for convection is similar to very 
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TABLE 1 

Summary of our detailed stellar evolution sequences. The 

COLUMNS give THE MODEL IDENTIFIER (S MEANS STERN, E IS 

EVOL, AND K IS KEPLER), the initial mass (Mq), the 

HELIUM CORE MASS PRIOR TO THE SECOND DREDGE-UP (Mq), THE 
HELIUM CORE MASS AFTER THE SECOND DREDGE-UP (Mq), 

INFORMATION THE END OF THE SIMULATION, AND THE FINAL FATE 
OF THE SEQUENCE ACCORDING TO OUR FIDUCIAL SAGB EVOL UTIO N 

PROPERTIES (mass LOSS, DREDGE-UP, AS DESCRIBED IN SeCT. I6.41| 
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fast mixing in semiconvective regions). 

The EVOL code has previously been used to study 
low-m ass (e.R. [ Hcrwig & Austin 2004) and massive AGB 
stars (iHerwid \2004mh ). KEPLER ha s in the past 
been applied to study massive stars (jWooslev et al.l 
l2002f ). but has not previously been used for AGB sim- 
ulatio ns. STERN has been used for low mass AGB 
stars (jLanger et al.lll999t iHerwig et al.ll2003t ISiess et alj 
2004f) as well a.s for massive stars ( Heger et al.l 120001 : 



Heger fc Langal[2000f ). 



3. PRE-AGB EVOLUTION AND THE INITIAL MASS 
RANGE FOR SAGB STARS 

In order to identify the processes that lead to SAGB 
star formation we calculate stellar evolution sequences 
with initial masses between 6.5 Mq and 13 M©, starting 
from the zero age main sequence until the completion of 
the second dredge-up or Ne ignition (Table ??). Up to 
the end of the second dredge-up, no mass loss is taken 
into account. The initial metallicity of our models is 
Z = 0.02. The effects of rotation or magnetic fields are 
not taken into account. 

3.1. H- and He-core burning 

The evolution o f stars toward the SAGB has been 
studied previously (Ritossa et al. 1996; Iben et al. 1997*; 
iGarcia-Berro et al.lll997 : Ritossa et al. 1999 : Siesa,,2006.) , 
and our simulations qualitatively confirm these results, 




Fig. 1. — Time evolution of convection zones and energy gener- 
ation for three evolution sequences with different mass, computed 
with STERN. The initial masses and evolution scenarios are: top 
panel: 5Mq, massive AGB, middle panel:11.5 Mq, SAGB, lower 
panel: 16.0 M© star, Fe-core, CCSN. Convective regions are indi- 
cated by green hatching while semiconvective regions are indicated 
by red cross-hatching. The energy generation from nuclear burning 
is shown in greyscale with a legend to the side in units of log erg 
g-i s-l. 

although quantitative differences occur. In our STERN 
models, a consequence of including semiconvection is 
that during core helium burning, a semiconvective layer 
limits the mixing between the inner helium burning core 
and the outer convective core, which still grows in mass 
(see also Fig. [1] below) . This decreases the lifetime of 
the core helium burning phase, because the available 
amount of helium is reduced, and leads to smaller he- 
lium and CO-core masses compared to models which use 
the Schwarzschil d crit erion for convection. 

Girardi et alj (|2000( ) studied the effect of convective 
overshooting on the maximum initial mass for which 
stars do not ignite carbon, M^p, and which defines the 
lower limit of SAGB stars. They find for models with- 
out overshooting a value of M^p of 6M0...7M0, while a 
moderate amount of overshooting reduces this by 1 Mq . 
In our models we find M^p = 7.5 Mq (EVOL/KEPLER), 
while our STERN models - without any overshooting - 
give Mup = 9.0 Mq. We will discuss these differences in 
the next paragraph. 

3.2. The second dredge-up 

The second dredge-up is a key differences between 
SAGB stars and massive stars that encounter Fe-core 
collapse. After core- He exhaustion, the core resumes con- 
traction while the envelope expands. As the star evolves 
up the asymptotic giant branch the envelope convection 
deepens, and eventually penetrates into the H-free core. 
Only due to this mixing event is the H-free core mass 
sufficiently reduced so that an electron-degenerate core 
can form which then cools and prevents Ne from ignit- 
ing. If the core mass after the 2nd dredge-up is smaller 
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than the Chandrasekhar mass, an electron-degenerate 
core will form and the He- and H-shells will eventually 
start the thermal pulse cycle. 

The dependence of the late evolutionary phases, in- 
cluding the second dredge-up, on the initial mass is illus- 
trated in the Kippenhahn-diagrams of three sequences 
computed with the STERN code shown in Figure [1] All 
models evolved through core-H and core-He burning. In 
the 5.0 Mq models, the hydrogen burning terminates, 
and the second dredge-up reduces the helium core mass 
by about 0.2 M0. When the helium shell source gets 
close to the bottom of the convective envelope, hydro- 
gen reignites, and the thermal pulse cycle starts. For the 
11.5 Mq model, central hydrogen and helium burning is 
followed by off-center carbon ignition. During the carbon 
burning phase the size of the helium core is reduced by 
a deep second dredge-up, after which the core becomes 
degenerate and thermal pulses develop. In the 16.0 M© 
case, convective core H- and He-burning is followed by 
core C-burning, and no 2nd dredge- up occurs. Ne ignites 
hydrostatically, and subsequent burning will lead to the 
formation of an iron core. 



In accord with previous work ( 


Ritossa et al.l Il996l: 


Iben et al.l 119971: iGarcia-Berro et al. 


19971: Ritossa et al.l 


19991: ISiessl'2006'). the second dredge- ud reduces the he- 



Hum core mass to values below the Chandrasekhar mass 
in our EVOL and STERN models. This leads to a clear 
definition of the upper mass limit of SAGE stars, as the 
critical mass between the occurrence and non-occurrence 
of the second dredge-up. For initial masses lower than 
this critical mass, Ne ignition is always avoided, while for 
initial masses higher than the critical mass the helium 
core mass is so large (~ 2.8 M©) that Ne always ignites. 
Our KEPLER models show a more complicated behav- 
ior: some show a second dredge-up depth which leaves 
helium cores with masses in between 1.4 M© and 2.8 M©. 
However, those with post-dredge-up helium cores above 
the Chandrasekhar limit all ignite core Ne burning. We 
conclude that a second dredge-up down to the Chan- 
drasekhar mass is required for an ECSN to occur, which 
thus defines our upper SAGB mass limit. This is also i n 
fine with the recent results of lEldridge fc ToutI ()2004aD . 
who do indeed find Ne shell flashes in some of their most 
massive models undergoing the second dredge-up; how- 
ever the dredge-up does proceed to the Chandrasekhar 
mass, and the suggested fate of these models is that of 
an ECSN. 

Figure [5] shows the helium core masses obtained in our 
detailed stellar evolution models. While the core mass 
after the 2nd dredge-up increases with initial mass for 
models computed with all three codes, differences arise 
with respect to the critical mass for second dredge-up. 
KEPLER and EVOL have similar final core masses, how- 
ever they differ slightly with respect to the maximum 
core masses. These differences are related to the treat- 
ment of convection and overshooting. STERN uses the 
Ledoux-criterion for determining the convective bound- 
aries, which naturally gives rise to smaller cores than the 
Schwarzschild-criterion. In these models no rotation was 
included, which - if included - would give significantly 
larger cores, due to rotationally induced m ixing during 
the hydrogen and helium burning phases (jHeger et al.l 
[2000.1 . 

Using EVOL and KEPLER we find the transition be- 
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Fig. 2. — Helium core masses for stars of various initial masses, as 
obtained using different stellar evolution codes (solid line: STERN, 
dashed line: KEPLER, dotted line: EVOL). The upper part of 
the line shows the maximum size of the helium core, prior to the 
second dredge-up. The lower part shows the size of the helium core 
just after the completion of the second dredge-up and prior to the 
onset of the TP-AGB. The light dashed horizontal line gives the 
lower limit for the final helium core mass for which the star may 
experience an electron-capture supernova. 

tween stars with and without a deep second dredge up 
at 9.25Mq. On the other hand, using STERN we 
find that stars more massive than 12 M© do not experi- 
ence a deep second dredge-up. The models in the mass 
range between 12 M© and 12.5 M© show a convective 
shell that develops on top of the helium burnin g layer 
(the so called dredge-out, c.f. IRitossa et al.lll999l ) which 
connects through a semiconvective layer with the bottom 
of the hydrogen-rich convective envelope. We find that 
protons are mixed into this hot layer and burn quickly. 
For a proper treatment of this interaction a scheme that 
solves the burning and mixing simultaneously is needed 
to follow the subsequent evolution of these stars. Since 
STERN is not equipped with such a scheme but cal- 
culates burning and mixing separately, we followed the 
subsequent evolution with very small timesteps until the 
code was not able to calculate further. Therefore, it 
remains unclear whether this semiconvective layer dis- 
solves and on what time scale. If it would, the helium 
core masses would be reduced to just below the Chan- 
drasekhar mass. If the semiconvective layer remains for 
the rest of the evolution of the star, it would allow Ne 
to ignite in the core and eventually lead to Fe-core col- 
lapse supernova. This renders the upper mass limit for 
SAGB stars according to the STERN models somewhat 
ambiguous in the range 12. ..12. 5 M©. 

In the EVOL models, the convective core overshoot- 
ing was calibrated by the observed width of the main 
sequence. However, stars rotate, and the STERN code 
usually takes this into account. The effect of rotation 
also widens the main-sequence and in this way STERN 
models with rotation ca n reproduce the observ ed main- 
sequence width as well (jHeger fc Langeill2000f ). In this 
study we compute non-rotating models with STERN, in 
order to avoid the complex question of how rotational 
mixing affects SAGB properties. As a drawback, the 
initial mass range for SAGB stars found with STERN 
is offset compared to the EVOL/KEPLER models by 
w 2.75 M©. Since the non-rotating STERN models do 
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Fig. 3. — Nuclear burning luminosity contributions and total 
luminosity as a function of time at the onset of the thermal pulse 
phase for a 9 Mq star calculated with STERN. At the time of 
second dredge-up hot-bottom burning starts, shown in the figure 
as a steep rise of the hydrogen luminosity. 

not reproduce the well-estabhshed main-sequence width, 
we prefer here the results of the EVOL /KEPLER models 
to derive the initial mass range for SAGB stars. Based 
on those, the upper mass limit for SAGB stars - and thus 
electron capture supernovae - is about 9.25 Mq (Fig- 
ure [2]). Without knowledge about the subsequent phase 
- the thermally pulsing SAGB phase - we can only say 
that the lower limit for electron capture supernovae is 
7.5 Mq, since stars with a lower mass do not ignite car- 
bon. 

4. THE TP-SAGB STELLAR EVOLUTION MODELS 
4.1. Thermal Pulses and Hot Bottom Burning 

Double shell burning of H and He on degenerate cores 
leads to periodic thermonuclear instabilities. These He- 
shell flashes or thermal pulses are an important site 
for nucleosynthesis in AGB stars, and cause mixing 
of the intershell region and — by way of the third 
dredge-up — mixing of processed inaterial to the surfac e 
(|Iben fc Renzinil[l98l IBusso et al.l[l999l: lHerwig|l200l . 
Thermal pulses of SAGB star s are similar to the rmal 
pulses of CO-core AGB stars f Ritossa et al.|[l996l ). In 
order to obtain quantitative information on these SAGB 
thermal pulse cycles, we calculate such model sequences 
for several initial masses (Table 1). 

As in massive AGB stars, most of the luminosity is 
produced by hot-bottom burning. During hot-bottom 
burning, hydrogen is transported convectively into the 
H-shell, and H-burning ashes are transported out of the 
shell into the envelope. In the more massive SAGB stars 
this hot-bottom burning starts immediately after the 
completion of the second dredge-up, and can proceed at 
very high temperatures. In our STERN models we ob- 
tain values of 1.0 x lO^K (10 Mq with = 1.16 Mq 
after 30 thermal pulses) and 1.1 x lO^K (11.5 M© with 
Mc = 1.27 Mq already at the first thermal pulse). The 
EVOL models show a similar trend with the 9.0 Mq 
model (E0099) reaching temperatures at the bottom of 
the convective envelope of 1.13 x lO^K after the 12th 
pulse. 

Hot-bottom burning could be stronger (or weaker) 
than in our calculations, e.g. due to convective overshoot- 
ing at the bottom of the convective envelope, or due to a 
larger (or smaller) convective efficiency than assumed in 



most MLT based stellar evolution calculations. In that 
case, the accretion of He on the core may be so much 
reduced that the core does not or only very slowly grow. 
We have performed some test calculations with enhanced 
convective extra mixing during the hot-bottom phase. 
These tests show so far a stationary H-shell, which slowly 
processes its envelope, and no core growth, which results 
probably in a massive ONe white dwarf. Whether this 
theoretical possibility is occuring in real stars is not clear 
because the physics of a convective boundary inside the 
H-shell is poorly known. 

SAGB stars show a He-peak luminosity during thermal 
pulses of around logL/LQ ~ 6, which is significantly 
lower than obtained in massive AGB stars which reach 
luminosities up to logL/LQ 8. This may explain why 
the third dredge-up is less efficient in terms of the dredge- 
up parameter A (see Sect. l4.2p . 

Extending the trend seen from low-mass to massive 
AGB stars, SAGB stars have smaller intershell masses 
(in the STERN 9 Mq model of 7 x 10"" Mq at a core 
mass of 1.06 Mq), and the interpulse time is also lower, 
ranging from 50 yr for the 11.5 Mq SAGB star with core 
mass of 1.27 Mq to 1000 yr for a SAGB star with a mass 
of 9.0Mq. 

4.2. Efficiency of the 3rd dredge-up 

The growth of the core during the TP-SAGB may be 
decreased by the dredge-up of material after a thermal 
pulse. The efficiency of the dredge-up (DUP) is expressed 
through the dredge-up parameter A = AMh/AA/dup, 
where AMjj is the core mass increase due to H-burning 
during the interpulse phase, and AMdup is the mass that 
is dredged up by the convective envelope. 

In the models calculated with STERN we did not ob- 
serve any dredge-up. This result is consiste nt with re- 
sults for non-rotating low mass AGB-stars (jSiess et al.l 
[200l from the same code which are also calculated us- 
ing th e Le doux-criterion for cori vection. iRitossa et al.l 
(|1996f) and ISiess fc Ptimol (I2006D find a sim ilar result. 
The recent models of iDohertv fc Lattanziol (120061 ) find 
very efficient dredge-up, e.g. A w 0.7 for a 9.5 Mq model. 
Observations clearly require a 3rd dredge-up in low-mass 
and massive AGB stars, since we see its result in terms 
of carbon and s-process enrichment in real AGB stars. 
However, the efficiency of the 3rd dredge-up in mas- 
sive and SAGB stars is not constrained observationally, 
partly probably because of the high dilution factor in the 
massive envelope. 

In order to get an idea about the efficiency of the 
3rd dredge-up in super AGB stars and the robustness 
of our and previous results, we studied the behavior 
of the thermal pulses also with the EVOL code. We 
calculated a 9Mq model (E0099) untfi the 12th pulse. 
This model was computed with a four times smaller fac- 
tor for the overshooting than the other EVOL models 
(/„^„. = 0.004) until the TP-AGB. This gives the star a 
smaller core than the regular models. On the TP-AGB 
a value of /over = 0.008 was used. The first thermal 
pulse starts after the completion of the second dredge- 
up, when the bottom of the convective envelope is at 
rrir — 1.17Mq. The surface luminosity after 12 pulses is 
logL/LQ = 5.07, the maximum helium luminosity dur- 
ing the 12th pulse logLnc/-^© — 6.17, and the duration 
of the interpulse period is ~ 500yr. 
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Fig. 4.— Evolution of the SAGB He-shell flash during the 12th 
pulse computed with the EVOL code (E0099). The dredge- up 
efficiency equals A = 0.5. 

After the eighth pulse, the ensuing mixing has the char- 
acteristics of a 'hot' dredge-up, first des c ribed f or massive 
low-metal licity AGB s t ars by iHerwid ()2004af ) and also 
found by IChieffi et all (pOOlf ) for Z = models. Any 
small amount of mixing of protons into the hot ^^C-rich 
layers — performed here by diffusive exponential over- 
shooting — leads to violent H-burning which increases 
the convective instability. Like a flame, this corrosive 
hydrogen burning enforces the penetration of the con- 
vective envelope into the intershell (see Fig. |4|) . For this 
situation, we find efficient dredge-up (A ~ 0.5), i.e. half 
of the interpulse core growth is dredged up, reducing 
the average pulse cycle core growth rate. Unfortunately 
the hot dredge-up phenomenon adds another source of 
uncertainty to dredge-up predictions as the dredge-up 
efficiency is extremely sensitive to the overshooting or 
extra-mixing efficiency at the bottom of the convective 
boundary. 

5. THE SAGB POPULATION SYNTHESIS MODEL 

Mass loss and the dredge-up are the two most impor- 
tant but also most uncertain processes that determine 
the final evolution of SAGB stars. Here we employ a 
simplified synthetic model that allows us to estimate the 
effect of different assumptions concerning these two pro- 
cesses on the initial mass range for ECSNe. 

5.1. A simple estimate 

We start with a simple back-of-the-envelope estimate: 
Stars that have, after carbon burning, a helium core 
mass larger than the Chandrasekhar mass (Mch) ex- 
plode as CCSN. The Chandr asekhar mass of a cold iron 
core is Mph-e M = 1.375 Mq (jSugimoto fc Nomotolll980l : 
lNomotol[l98ll ). In order to form an ECSN, the core mass 
has to grow from the mass at the beginning of the TP- 
AGB, Mc{2DUP), to the Chandrasekhar mass by 



AMe = Afch - AU2DUP). 



(1) 



This value depends strongly on the initial mass as Fig- 
ure [5] shows. 

Whether the core is able to grow by this amount de- 
pends only on the mass of the envelope, the core growth 
rate and the mass loss rate. Given these quantities, 
AMc^max is the maximum mass that the core can grow. 

The core growth rate due to nuclear burning is dAfc/dt. 
Due to the 3rd dredge-up, the value for the core growth 
rate can decrease. We correct for this by introducing a 
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Fig. 5. — Mass AA/c by which the core needs to grow during the 
TP- AGB in order to reach the Chandrasekhar mass, as function 
of the initial mass. The upper boundary is given by stars that 
do not have a second dredge-up, so their cores are larger than the 
Chandrasekhar mass. 

factor 1 — A. The timescale on which the envelope of the 
star will be lost is 



dM/dtc 



(2) 



and multiplied by the core growth rate this gives the 
maximum value that the core can grow. This gives an 
approximate relation for the growth of the core during 
the TP-SAGB phase 



dAf/dtc 



X (1 - A)dA/c/dt. 



(3) 



For a typical, but constant, core growth rate of Mc = 
5 • 10~^ M0yr~^ , and an envelope mass of A^cnv = 10 M©, 
Fig.[n]shows AA^cmax as a function of the mass loss rate, 
for two different values of A (no dredge-up and A = 0.9). 

Figure [5] shows that in order to have an initial mass 
range for ECSN of, for example, IM©, the core growth 
during the SAGB phase must be of the order 0.1-0.2 M©. 
Figure [6] shows that if SAGB mass loss is larger than 
« lO~^M0/yr such a core growth can not be achieved, 
even for inefficient 3rd dredge-up. For mass loss rates 
below ~ 10~^MQ/yr, however, a core growth of a few 
0.1 Mq is predicted even if A = 0.9. Compared with 
the em pirical mass loss rates derived bv Ivan Loon et aTl 
(|2005l hereafter L05, c.f. Tablell]) it is clear that the intial 
mass range for ECSN is sensitive to the third dredge-up. 

5.2. Synthetic SAGB evolution 

A quantitative estimate of the initial mass range for 
ECSN can be obtained through a s ynthetic model for 
the TP- AGB phase, similar to that of llzzard et"all (|2004 
hereafter 104) for AGB stars, which is based on d etailed 
AGB models from Karakas (jKarakas et al.l 120021 ). The 
extension to SAGB stars is made by fitting the TP- AGB 
evolution of detailed stellar evolution models (STERN) 
presented above, specifically over the mass range between 
7 and 11.5 M© in initial mass. 

Based on the SAGB STERN evolution sequences with 
up to 30 thermal pulses, we derive fits for luminosity, 
radius, and Q-factor (see Sect. 15.23]) . as function of the 
core mass (A/c), the envelope mass (Afonv) and as sec- 
ondary parameters the metallicity {Z) and the envelope 
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Fig. 6. — Expected core growth during the TP-SAGB (Eq.[3} as 
function of the mass loss rate, for two different dredge-up efficien- 
cies A as labeled. A constant mass loss rate, an envelope mass of 
10 Mq, and a constant core growth rate of 5 X lO~^M0/yr~l are 

Ji§FtiH:#!n abundance (Xh). Since the SAGB evolution 
models have entered into a quasi-steady state regime, 
these fits are good approximations for the subsequent 
evolution of SAGB stars during the TP-AGB in mass 
(total, core and envelope), luminosity and radius. 

We then use these analytic expression as basis for 
our synthetic TP-SAGB model. As starting values for 
our synthetic SAGB calculation we use total mass, core 
mass, and envelope hydrogen abundance after the sec- 
ond dredge-up. First the luminosity is calculated from 
the initial core and envelope mass, then the radius is cal- 
culated, which is a function of the previously calculated 
luminosity and the envelope mass, and finally the core 
growth is calculated and integrated over a timestep dt. 
From these quantities, the effective temperature, mass 
loss rate, the resulting new mass of the envelope, and 
the new mass of the core are calculated. The new core 
mass and envelope mass are used as input for the next 
timestep. 

In following subsections we describe the b asic outline of 
our s ynthetic model (for details we refer to llzzard et al.l 
[2001 . 

5.2.1. Luminosity and Radius 

We follow 104 and model the luminosity with two 
terms, one which contains a core-mass-luminosity rela- 
tion (CMLR) and one term due to hot-bottom burning. 
The total luminosity of the star can now be written as 
(cf. their Eq. 29) 

L = fdiftLcMLR + Lcnv) L©, (4) 

where Lcmlr is the core mass-luminosity relation given 

by 

LCMLR = 3.7311 10* X 

max[(Mc/ Mq - 0.52629) (2.7812 - M^/ Mq), 

1.2(Me/M0-O.48)] (5) 

if the core mass at the first thermal pulse, Ale itp is 
>0.58Mq. 

Lenv is the contribution due to hot-bottom burning 
(e.g., 104:32), 

1.3 r 



Lpnv= 1.5010^ 



1 + 0.75 1- 



Z 
(102 



X max 



Mc 
Mo 



1 AM 



c.nodup 



M, 



- 0.75 



o 



(6) 



Mc is the core mass, M^nv is the envelope mass. 
AMc,nodup is the change in core mass without third 
dredge-up and is defined by AMc,nodup = Mc,nodup - 
Afc,iTP with Mc_nodup the core mass as if there was no 
third dredge-up and Mc.itp the core mass at the first 
thermal pulse. Z is the metallicity. Note that we use a 
lower exponent than 104 in the contribution of the enve- 
lope mass, i.e. 1.3 instead of 2, which resulted in good 
fits for models between 7 and 11.5 M©. 
The function 



ft = min 



J^AMc^„odup/Mo 



0.04 



0.2 



,1.0 



(7) 



accounts for the steep rise in luminosity at the beginning 
of the TP-AGB. The function 

/d = 1 - 0.2180exp [-11.613(Afc/ M© - 0.56189)] (8) 

corrects for the short timescale dips in the luminosity 
during the thermal pulse cycle. 

For the fit to the radius we use an expression of the 
same form as given by 104, but with coefficients adjusted 
to the STERN models: 



log(/ri?) 
with 



-0.26 + 0.75 log (L/ Lq) - 0.41 log Afo„v/ M8) 



A = 0.091og(Menv/Afc„v,0) 



(10) 



a factor that accounts for the removal of the envelope, 
where Afcnv.o is the mass of the envelope at the first 
thermal pulse. This correction factor is determined by a 
fit to a 9 Mq model to which an extreme mass loss rate 
of lO^'^Mo/yr was applied. For envelope masses below 
Mcnv = 2 Mq the fit predicts too large radii and it is not 
valid for temperatures below 2500 K as calculated from 
R and L. 

5.2.2. Third Dredge-up 

For the dependence of the thir d dredge-up on t he in i- 
tial mass we use the data from iKarakas et al.l ([2002) . 
Our own EVOL SAGB models, however, show smaller 
dredge-up (A = 0.5 for Min i = 9M q) than the extrapo- 
lation of the IKarakas et all (|2002f ) data (Figure O. We 
therefore extend the fit to higher masses with a relation 
that reflects our own data at Mini = 9 Mq . To simulate 
a situation with no dredge-up, we also include in our 
synthetic code an option to set A = 0. 

5.2.3. Core growth and hot-bottom burning 

The growth rate of the He core in the inter-pulse phase 
is given by 



dMe 
dt 



X L 



(11) 



where L is the total luminosity of the star, and Q gives 
the mass of nuclear ashes accreted onto the core per en- 
ergy released by the star. Q depends on several model 
properties, especially on the hot-bottom burning effi- 
ciency, and less on the chemical composition of the en- 
velope. For massive AGB and SAGB stars its strength 
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Fig. 7. — Dredge-up efficiency as a function of init ial mass. The 
dotted line extrapolates the data by Kar akas et al.l ((20021) , while 
the solid line gives the modification based on our SAGB dredge-up 
stellar evolution sequence. 




Fig. 8. — Ashes accreted on core per unit stellar energy release 
(Q) as a function of the mass of the envelope. For envelopes less 
massive than of 4 Mq there is no hot-bottom burning. 

depends on the envelope mass. If the hot-bottom burn- 
ing is efheient, Q can be smah (see discussion in § I4.ip . 
and the core may not significantly grow at all. Fig. [5] 
shows the decrease of Q with increasing envelope mass. 
We parameterize Q as 



g = min[1.43 x 10"", 1.40 x 10"" 



4.166 X 10 



-12 



1.5 X 10 



-12 -^^cnv 1 
Mr. ^ 



(12) 



This parameterization is in reasonable ag reement with 
104 wh o set Q to 1.585 x 10"" M© L^Vr- MH^rlev et all 
(I2OOOI) who found 1.27 X 10"" M© L-Vr"\ and 1.02 x 
10~^^ Mq Lg^yr"^ in iWagenhuber fc Groenewegenl 



5.2.4. Mass loss 

As discussed in § 15.11 the mass loss of SAGB stars is 
the other important factor to determine the initial mass 
range of ECSN. SAGB stars are 0-rich (because of hot- 
bottom burning), and have stellar parameters around 
log Toff = 3.5 and logi/L© = 5 at solar metallicity. It 
is not clear what the dominant mass loss mechanism for 



TABLE 2 

Mass loss rates for SAGB stars with a typical value for 

THE luminosity OF logL/ Lq = 5 AND AN EFFECTIVE 

temperature of T(,ff = 3000K. 







Type 






Rate 


Reimers [t] 


= 1) 


Red Giants 


~ 5 


10" 


"Me/yr 


Reimers {rj 


= 4) 


Red Giants 


~ 2 


10" 


SMe/yr 


Schroder & 


Cuntz 


Super Giants 


1 


10- 


SMe/yr 


van Loon 




AGB/RSG 


- 3 


10- 


SMe/yr 


Blocker 




AGB 


- 6 


10- 


^Me/yr 


Vassiliadis 


Wood 


AGB 


4 


10- 


SMe/yr 



these stars is. Are they cool enough to develop dust- 
driven winds or is mass loss simply driven by radiation 
pressure? 

Table [2] shows a compilation of observational and theo- 
retical mass loss rates, for a combination of typical SAGB 
parameters. Note, however, that these rates are not con- 
stant over time, and that the variation itself during the 
AGB phase is important for the final outcome. We pref- 
erentially use the observed mass loss rates for massive 
AGB stars and red supergiants by Ivan Loon et all (|2005l 
L05). If dust- formation does not pl ay an importan t role, 
then the Reimers mass loss rate (jReimer may 
be applicable. It is derived from observations of RGB 
stars with a small range in tem peratures and radii, how- 
ever, ISchroder fc Cunta (|2005| ) have revised the Reimers 
rate. For the more massive RSG stars their new ap- 
proach, which also includes surface gravity, gives about 
three times larger mass loss than the Reimers formula. 
This places it, for given temperature and luminosity (c.f. 
Tabled]), within a factor of 2 of the observational mass 
loss det ermina tion b y L05 . The mass loss formula by 
IVassiliadis fc WoodI (Il993l) (hereafter VW93), which is 
often used for AGB star evolution calculations, is also 
close to the o bservational va lue. The AGB mass loss 
formulated by Blo cker! (jl995l), and bas ed on the hydro- 
dynamic wind models bv iBowen! (1988) has a higher lu- 
minosity exponent, and gives very high mass loss rates 
for SAGB stars. From our first estimate in Sect. IS.ll it is 
clear that with the Bloecker mass loss SAGB stars would 
never explode as ECSN. 

6. RESULTS 

We perform a series of synthetic calculations, with two 
assumptions on third dredge-up and three assumptions 
on mass loss. For dredge-up we assume either the pa- 
rameterization of § 15.2.21 or A = 0. For mass loss, we 
consider the cases Reimers, L05, and VW93 f§ I5.2.4p . 

6.1. Initial mass range for ECSN 

The resulting initial mass ranges for ECSN are illus- 
trated in Figure[9]for the case with parameterized A, and 
in Figure [TUl for A = 0. Stars that end their evolution as 
white dwarf, i.e. below the Chandrasekhar mass, do not 
explode as ECSNe. With the parameterized prescription 
for the third dredge-up, the width of the initial mass 
window for which ECSN occurs is between 0.25MQand 
0.65 Mq, depending on the assumed mass loss rate. The 
mass loss prescriptions of L05 and of VW93 give an initial 
mass window of 0.20 — 0.25 M©. With zero dredge-up, 
the core grows at the maximum possible rate. The width 
of the initial mass window for ECSN is between 1.4 M© 
for the Reimers mass loss rate and 0.45 — 0.5 M© for the 
VW93 and L05 mass loss rates. 
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TABLE 3 

Initial mass limits for ECSN and ratio of ECSN to SN as a function of the dredge-up efficiency and mass loss prescription. 





A = 


parameterized 






A = 






Afiow/Me 




% EC 


Miow/M0 


Mhish/ Mq 


% EC 


Reimers (rj = 4) 


8.67 


9.25 


8.4 


7.86 


9.25 


19.7 


VW93 


9.03 


9.25 


3.2 


8.82 


9.25 


6.2 


LOS 


9.00 


9.25 


3.6 
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7.1 
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Fig. 9. — Final core mass as a function of initial mass based 
on synthetic SAGB calculations. The solid line indicates the post 
second dredge-up core mass, the short dashed line indicates the 
final core mass using the Reimers mass loss rate (rj = A), the dashed 
line using the L05 mass loss rate, and the dash-dotted line using 
the VW93 mass loss rate. The shaded region indicates the initial 
mass range for ECSNe for the LOS mass loss rate. 



1.45 



1.4 



S 1.35 



1.3 



1.25 



1.2 



■ post 2DUP 

- Reimers 


' ' 1 ' ' 




'_ van Loon 






■ Vassiliadis & Wood 












/ / 
/,' 

/ 

/' / 
/ 






/ 






// / 
// / 
/,' / 






/ 

. A . . \ . . /. . \ . . 


. . 1 , , 





7.5 



8 8.5 



9.5 



Fig. 10. — Same as Figure |9] but for calculations assuming no 
dredge-up. Note different x-axis scale. 



6.2. ECSN fraction 

Based on the inferred mass ranges from the synthetic 
model, we determine the ratio of the number of ECSNe 
to the total number of SNe. Table [3] gives an overview of 
the results for the cases of parameterized dredge-up and 
without dredge-up (A = 0), assuming the Salpeter IMF. 
The value of A has a strong influence on the predicted 
fractions. With the parameterized dredge-up and the 
VW93 or LOS mass loss rates the ECSN fraction of all 
supernovae is about 3.5%. With the Reimers mass loss 
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Fig. 11. — Final mass of the remnant as a function of the initial 
mass. Remnant regimes are shaded as light grey: CO-white dwarf; 
grey: ONe white dwarf; white: ECSN; dark grey: CCSN. The final 
mass is either the WD mass or the stellar mass at the time of SN 
explosion. The dashed line indicates the line of initial mass equal 
final mass. 

rate 8% of all supernovae are ECSN. The largest ECSN 
fraction of 20% is obtained without dredge-up and using 
the Reimers mass loss rate. 

6.3. Final masss and SN properties 

Figure [TT] shows the initial-final mass relation for the 
mass range from 1.0 to 14 using the parameter- 
ized dredge-up prescription and the LOS mass loss rate. 
For ECSNe, we find a large spread in progenitor and 
envelope masses. The least massive SAGB SN progeni- 
tors lose almost their entire envelope, growing the core 
just barely enough to still make an electron capture su- 
pernovae before the envelope is lost. The most massive 
SAGB SN progenitors, on the other hand, undergo very 
little TP-AGB mass loss before they explode, and con- 
tain a massive hydrogen-rich envelope at that time. 

This diversity is a natural consequence of the com- 
petition between core growth and mass loss during the 
SAGB stage, and thus independant of the choice of 
mass loss rate and dredge-up parametrisation. The ex- 
pected envelope mass range, from almost zero to about 
8 Mq (Fig. [12]), implies a diversity of supernova hght 
curves of ECSNe, which may range from light curves of 
so called Type lib and Type IIL supernovae to those of 
typic al Type IIP supernovae (jFalk fc Arnettlll977l : lYounel 
I2OOI . 

However, ECSNe may show three properties which 
might allow to distinguish them from ordinary Type II 
supernov ae. First, they might produce low-energy ex- 
plosions (Kitaura et al. 200( 1) and p ossibly low neutron 
star kicks ( Podsiadlowski et al.ll200-^ . A consequence of 
the low explosion energy may be a small nickel mass pro- 
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Fig. 12. — Mass of the envelope just prior to the explosion as 
function of the initial mass for two mass loss rates. 
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Fig. 14. — Mass loss rates according to LOS as function of initial 
mass, for the models shown in Fig |13l 
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Fig. 13. — Maximum luminosity as funtion of the initial mass at 
the end of the evolution. Progenitors of CO white dwarfs reach 
luminosities up to logL/L© ~ 4.6, progenitors of ONe white 
dwarfs reach luminosities up to logL/LQ ~ 5. Progenitors of 
ECSN are the most luminous, with log L/Lq > 5, while progeni- 
tors of the least massive CCSN have pre-explosion luminosities of 
log L/ Lq ~ 4.6 



duced by the explosion, and thus a low luminosity of the 
tail of the light c urve which is produc ed by the decay 
of ^^Ni and s^Co (|Kitaura et all [2006^. The Type IIP 
SN 1 997D may provide an example ( Chuga i fc Utrobiiil 
|2000| ). Second, however, the enormous mass loss rate 
of the supernova progenitor (Fig. star may produce 
clear signatures of a supernova-circumstellar medium in- 
teraction in the supernova light. Such signatures are 
in particula r exceptionally a bri ght and long-lasting 
light curve (jSoUerman et all |2001[ ) , and narrow hydro- 
gen e mission lines supe rimposed to a typical SN II spec- 
trum (jPastorello et all [20021. Third, ECSN progenitors 
are extremely bright, with luminosities of the order of 
10^ L0 (Fig. [121). Thus, p rogenitor i dentifi c ations on pre- 
explosion images (|Maund fc Smarttll2005t iHendrv et all 
[2006) might be able to identify ECSNe. They may be dis- 
tinguished from very massive (> 20 M©) progenitors of 
similar luminosity by their much cooler effective temper- 
atures (< 3000K for ECSN progenitors versus ~ 3400K 
for CCSN progenitors). 
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1x10' 2x10 3x10 

time on TP-AGB / yr 
15. — Synthetic AGB mass evolution for an initial mass of 
Upper panel: mass loss rate for two cases (LOS and VW93, 
15. 2.41 1: the superwind regime in the VW93 mass loss prescription 
sets in at lOOOyr right at the beginning of the AGB phase. Lower 
panel: evolution of the core and the total stellar mass. 

6.4. The reference model, examples 

As shown above, the results of our synthetic SAGB 
calculations employing the VW93 and the LOS mass 
loss prescriptions are rather similar. Since it is unclear 
whether the Reimers mass loss rate is really applicable, 
and as the L05 mass loss rate relies on very recent ob- 
servations, we adopt the latter as the fiducial mass loss 
prescription for our synthetic SAGB modeling. Concern- 
ing the third dredge-up efficiency, we adopt the mass 
dependant formulation shown in Fig. [7| as our reference 
efficiency. These assumptions define our reference model 
for the synthetic SAGB evolution. 

In the following, we discuss some explicit examples to 
illustrate the TP-SAGB evolution, and to further moti- 
vate the choice of our reference model and analyse its 
uncertainty. The first example is a star with an initial 
mass of 9.1 Mq, with a He-core mass at the end of the 
second dredge- up of 1.348 M©. 

During the evolution on the TP-SAGB the luminosity 
first increases from logi/L© — 4.9 to logL/L© = 5.02, 
and then drops slightly due to decreasing envelope mass 
which decreases the efficiency of hot-bottom burning. As 
a result, the inter-pulse core growth rate increases from 
3 X lO-'^Me/yr to 1.5 x lO^'^Mg/yr, but the effective 
core growth is significant lower (by about factor 0.5) due 
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time on TP-AGB / yr 
Fig. 16. — Similar to Figure [Hbut for a 8.8 Mq star. The core 
does not reach the Chandrasekhar mass but loses its envelope and 
becomes a massive ONe white dwarf. 

to the effect of dredge-up. The mass loss rate increases 
from 3 X lO^^Me/yr to 1 x lO^^Mo/yr (Fig. [151). In 
this model the SAGB ends after ^ 4.4 x 10 yr when 
the core reaches its Chandrasekhar mass. The remaining 
envelope has a mass of 4.27 Mq. 

For a 10% larger dredge-up efficiency the SAGB time 
increases to ^ 4.8 x lO'^yr, and the remaining envelope 
decreases to 3.82 M©. For a 10% smaller dredge-up ef- 
ficiency the SAGB decreases to ~ 4.0 x lO^yr, and the 
final envelope mass increases to 4.61 Mq. For the case of 
no dredge-up the SAGB time is ~ 2.4 x IC^yr and the 
final envelope mass is 6.1 Mq. 

Using the mass loss prescription of VW93 and the pa- 
rameterized dredge-up, the result is only different in the 
mass loss history, resulting in a smaller final mass. The 
main reason is that this mass loss prescription accounts 
for the superwind phase for TP-AGB stars, which gives 
a significantly different evolution of the envelope. The 
TP-AGB phase starts with low mass loss but makes a 
transition to the superwind phase after 1000 yr with mass 
loss rates around 1 x 10^** Mq/ yr, decreasing slowly due 
to the waning luminosity (Fig. [151 Upper Panel). Af- 
ter 4.7 X lO'^yr, the core reaches its Chandrasekhar mass 
with a final envelope mass of 2.43 Mq. 

An initially 8.8 Mq star with the same mass loss and 
dredge-up (Fig. [TBI) becomes an ONe WD. The He- 
core mass at the beginning of the SAGB is 1.296 Mq. 
During the SAGB evolution the luminosity increases to 
logL/Lg = 4.95, then decreases due to mass loss. As- 
suming the LOS mass loss rate, the envelope is lost after 
'--^ 1 X 10^ yr. The mass loss rate increases steadily dur- 
ing the TP-AGB phase, but as the star reaches a surface 
temperature of 2500K (Sect. 15. 2p we assume a constant 
mass loss rate of M = 1 x 10^''MQ/yr during the last 
phase of the evolution. This is consistent with the ob- 
servations of L05 who find only one star with a higher 
mass loss rate. The remaining ONe core has a mass of 
1.338 Mq. Repeating this calculation with the mass loss 
rate of VW93 does not show large differences (Fig. [TBt . 

7. CONCLUDING REMARKS 

We show that both the lower initial mass for SAGB 
stars and the maximum initial mass for ECSNe sensi- 
tively depend on the assumptions for mixing during core 
H- and He-burning. EVOL models which include core 
overshooting, and which are consistent with the observed 
width of the main sequence, predict a smaller maximum 



initial mass for SAGB stars. Rotation would act simi- 
lar to overshooting during the core burning phases. The 
STERN models include neither rotation nor overshoot- 
ing, and the maximum initial mass for SAGB and EC- 
SNe is larger by up to 2.5 Mq. Equally important is the 
treatment of semiconvection during He-core burning. 

On the other hand, the lower initial mass limit for EC- 
SNe is determined by the stellar properties on the SAGB. 
Most important are the third dredge-up efficiency, the 
mass loss rate, and the hot-bottom burning efficiency 
and its dependence on the adopted convection theory for 
the envelope. In general, larger mass loss, larger dredge- 
up efficiency and large hot-bottom efficiency all decrease 
the initial mass range for ECSNe, or even suppress the 
ECSN channel. In order to increase the accuracy of the 
transition initial mass between ECSNe and CCSNe and 
of the lower mass limit for ECSNe, these classical issues 
of stellar evolution need to be improved specifically for 
the initial mass range of 6 to 12 Mq. 

We have discussed here the SAGB stars with C-ignition 
and formation of ONe cores as the most likely progenitors 
of an ECSN class of supernova. It is in principle possible 
that initially less massive stars that develop CO cores 
could increase their core size on the TP-AGB up to the 
Chandrasekhar mass resulting in a supernova explosion. 
Despite the uncertainties still involved we can rule out 
this possibility for solar metallicity. There are two main 
reasons that prevent these SN1.5a from occurring: First, 
the mass loss would have to be much lower than observed. 
Second, models predict that the third dredge-up is larger 
for massive AGB stars with initial mass between 4 an d 
7Mo than for SAGB stars rFig.[7l iKarakas et al.|[200l . 
This makes it even more unlikely for massive AGB stars 
to significantly grow their cores. 

We did not take into account mass loss until the be- 
ginning of the thermally pulsing phase. If mass loss were 
applied during the main sequence and up to the TP- 
AGB, less than half a solar mass would have been lost 
(jSiesj 120061 their Table 5). This may shift the quoted 
initial masses to a slightly higher value. 

We note that there is large disagreement between the 
differe nt studies on the d r edge- up efficiency in SAGB 
stars (iRitossa et all 119961 : iDohertv fc Lattanziol l2006l : 

ISiess fc Pumol I2006D . Whereas IRitossa et al.l (|1996[ ) 

and iSiess & : Pumo' ('20 06f) find negligi ble amounts of 
dredge-up, Doherty & Lattanziol (|2006( ) find very effi- 
cient dredge-up with A ~ 0.7. This situation is sim- 
ilar to the divergent modeling results obtained in the 
past on the third dredge-up in low-mass AGB stars. For 
the low-mass regime there is now some consensus that 
most of the differences where related to different phys- 
ical and numerical treatments of the convective bound- 
aries. Possibly the same applies to the divergent results 
for SAGB stars. In our synthetic model we adopted a pa- 
rameterized prescription (Fig.[7l) that is based on state- 
of-the-art full stellar evolution calculations. To test the 
effect of dredge-up we also considered a case with no 
dredge-up. Clearly, the third dredge-up in SAGB stars 
needs to be studied in more detail. It is closely related 
to the mixing conditions at the bottom of the convec- 
tive boundary. This is a hydrodynamic situation which 
requires multi-dimensional simulation which is compli- 
cated by the fact that for these extremely mas sive cores 
the dredge-up seems to be hot (jHerwi gl [2004al ). I.e., any 
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TABLE 4 

Final core and envelope masses for different initial masses. The first three columns give the initial conditions of the 
MODELS. Models marked with an asterisk are calculated until the end of the second dredge-up with EVOL (Table ??), the 
other models are interpolated assuming a linear relation between the initial mass and the core mass at the end of the 
second dredge -up. T he fourth and the fifth column give the final c ore an d envelope mass for two different mass loss 

prescriptions 15. 2. 411 AND FOR A PARAMETERIZED DREDGE-UP EFFICIENCY 15. 2. 211 . THE SIXTH AND SEVENTH COLUMN GIVE THE FINAL 
CORE AND ENVELOPE MASS FOR A DREDGE-UP EFFICIENCY OF ZERO. In CASE A COLUMN CONTAINS "WD" OR "NS", THE REMNANT IS 
RESPECTIVELY' A WHITE DWARF OR A NEUTRON STAR. AlL MASSES ARE IN UNITS OF A SOLAR MASS. 



initial conditions 


final conditions param 


A 


final conditions A = 











LOS 


VW 93 


LOS 


VW 93 


A/tot 


Afc 




A/c 




A/c 


A/onv 


A/c 




Afc 


A/cnv 


6.5000* 


0.9900 


5.5100 


1.0039 


WD 


1.0060 


WD 


1.0816 


WD 


1.0887 


WD 


7.5000* 


1.0700 


6.4300 


1.0886 


WD 


1.0903 


WD 


1.1734 


WD 


1.1739 


WD 


8.5000* 


1.2440 


7.2560 


1.2809 


WD 


1.2761 


WD 


1.3358 


WD 


1.3255 


WD 


8.7500 


1.2875 


7.4625 


1.3282 


WD 


1.3223 


WD 


NS 


0.2075 


1.3654 


WD 


8.8000 


1.2962 


7.5038 


1.3376 


WD 


1.3316 


WD 


NS 


1.3774 


1.3736 


WD 


8.9000 


1.3136 


7.5864 


1.3563 


WD 


1.3504 


WD 


NS 


3.2631 


NS 


1.4514 


9.0000 


1.3310 


7.6690 


NS 


0.0221 


1.3693 


WD 


NS 


4.7723 


NS 


3.0680 


9.1000 


1.3484 


7.7516 


NS 


4.2672 


NS 


2.4311 


NS 


6.0628 


NS 


4.7208 


9.2000 


1.3658 


7.8342 


NS 


6.7800 


NS 


5.7643 


NS 


7.2711 


NS 


6.6314 


9.2500 


1.3745 


7.8755 


NS 


7.8321 


NS 


7.7731 


NS 


7.8481 


NS 


7.8138 



small amount of H that could be mixed across the con- 
vective boundary will instantly burn violently with all 
the associated feedback on the evolution of the convec- 
tive instability in that region. On the other hand, the 
amount of material that is dredged up from the He in- 
tershell will be very small, even with large values of A, 
due to the thin intershell, and will therefore easily dilute 
in the envelope. 

It is presently not known whether ECSNe from 
SAGB stars contribute to the r-process pattern in the 
universe. Explosions from stars in this mass range have 
been investigate d as promi sing site for the astr ophys- 
ical r -process (iWheeler et al.l Il998t ISumiyoshi et alJ 
20011: IWanaio et all l2003l) . based on the work of 



Hillebrandt et al.l (I1984D who ex ploded a ONe 
model from iNomotol (|19M [l98l . Other groups were 
not able to confir m the r-process contribution due to the 
lowe n tropv fe.g..lBurrows fc Lat timer 1985: B aron et al.l 
19871: lMa\de fc WilsonI 119881 : iBethe & WilscmrfT985i r 
Kitaura et al.l (|2006l ) ruled out this possibility based on 
updated physics and two different nuclear equations of 
state. 



In any case, our study outlines that ECSNe from SAGB 
stars are likely to occur, if only at a level of a few per- 
cent of the Type II supernova rate in the local universe. 
However, at low metallicity, the key physical ingredients 
to the evolution of thermally pulsing SAGB stars may 
change. In particular, the stellar wind mass loss rate 
may be lower, which might open the ECSN channel ap- 
preciably, and may even allow Type 1.5 supernovae. This 
issue will be discussed in a forthcoming paper. 
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